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Abstract. We have computed the dust temperature distribution to be expected in a pre-protostellar core in the 

^ ' phase prior to the onset of gravitational instability. We have done this in the approximation that the heating 

1/^ , of the dust grains is solely due to the attenuated external radiation field and that the core is optically thin to 

•/< ■ its own radiation. This permits us to consider non spherically symmetric geometries. We predict the intensity 

>— ' ' distributions of our model cores at millimeter and sub-millimeter wavelengths and compare with observations 

of the well studied object L1544. We have also developed an analytical approximation for the temperature at 

the center of spherically symmetric cores and we compare this with the numerical calculations. Our results show 

(in agreement with Evans et al. 2001) that the temperatures in the nuclei of cores of high visual extinction (> 

■~^ ■ 30 magnitudes) are reduced to values of below ~ 8 K or roughly half of the surface temperature. This has the 

i-i^ ' consequence that maps at wavelengths shortward of 1.3 mm see predominantly the low density exterior of pre- 

Mh^ protostellar cores. It is extremely difficult to deduce the true density distribution from such maps alone. We have 

Q ■ computed the intensity distribution expected on the basis of the models of Ciolek & Basu (2000) and compared 

5— ( ' with the observations of L1544. The agreement is good with a preference for higher inclinations (37° instead of 

f/3 , 16° ) than that adopted by Ciolek & Basu (2000). We find that a simple extension of the analytic approximation 

C^ ' allows a reasonably accurate calculation of the dust temperature as a function of radius in cores with density 

distributions approximating those expected for Bonnor-Ebert spheres and suggest that this may be a useful tool 

for future calculations of the gas temperature in such cores. 

X, 

M Key words. Molecular Clouds - Dust Grains 

1. Introduction to derive the density distribution in such sources from ob- 

servables such as the millimeter emission and the mid-IR 



The structure of the high density cores embedded in , ,. / ^^ , , nnnn\ 

. / . absorption (e.g. aacmann et al. 2000) 

nearby molecular clouds is currently a matter oi great in- 
terest. The fate of a protostar depends sensitively on the 
initial conditions prior to collapse and it is difBcult to infer 
this from studies of "protostars" where a luminous object 
has already formed. Considerable observational effort has 
therefore been expended in the attempt to detect a "pre- 
protostellar core" , which is taken to be a cold dense core 
whose column density and internal pressure exceed consid- 
erably those found in the surroundings. The high density 
nucleus of the dark cloud L1544 is an example of such an 
object (Tafalla et al. 1998, Ward-Thompson et al. 1999, 
WiUiams et al. 1999, Ciolek & Basu 2000a,b) and in fact 
the models of Ciolek & Basu (2000a, CB hereafter) sug- 
gest that it will become unstable in around 30000 years 
from now. However, these inferences depend on our ability 



In particular, the emission of dust grains at millime- 
ter and submillimcter wavelengths is sensitively dependent 
not only on the grain column density but also upon the 
grain temperature. In clouds such as L1544, the grain tem- 
perature is known to be of order 12 K on the basis of the 
observed spectral energy distribution (Andre et al. 2000) 
and it has been normal practise to assume such regions 
to be isothermal. However, it has been known for some 
time that one can expect appreciable temperature gra- 
dients within cores heated by the (external) interstellar 
radiation field (e.g. Leung 1975; Mathis et al. 1983, here- 
after MMP) as well as temperature differences between 
grains of differing optical properties. For example, for sil- 
icate grains, MMP predict grain temperatures as low as 
6 K in the nucleus of cores of visual extinction ~ 50 mag- 
Send offprint requests to: M. Walmsley nitudes. Such low temperature grains will contribute neg- 
Correspondence to: walmsley@arcetri.astro.it ligibly to the emission at wavelengths below 1 mm and 



A. Zucconi et al.; The dust temperature distribution in prestellar cores 



this suggests that a new study of grain temperatures in 
dense cloud cores similar to LI 544 is warranted. 

Another characteristic of the observed pre-protostellar 
cores is that they clearly show large departures from 
spherical symmetry in both the maps of dust emission 
and absorption (Bacmann et al. 2000). This is often in- 
terpreted as being due to flattening along the direction of 
the mean magnetic field resulting in a disk-like configu- 
ration (Ciolek & Basu 2000a, b). Furthermore, as shown 
recently by Galli et al. (2001), a cloud core modeled as a 
thin disk perpendicular to the large-scale magnetic field 
and supported against its self-gravity by a combination 
of gas pressure and magnetic forces, is not necessarily ax- 
isymmetric. Direct evidence for the importance of a mag- 
netic field in determining the structure of L1544 comes 
both from polarization measurements at 850 fiia (Ward- 
Thompson et al. 2000) and observation of OH Zeeman 
splitting (Crutcher & Troland 2000). The magnetic field 
inferred from the Zeeman measurement is consistent with 
the model predictions of CB but this result is sensitive 
to the assumed inclination of the magnetic field relative 
to the plane of the sky (16°, according to CB). It is rele- 
vant also that the observed direction of polarization is not 
consistent with naive CB model expectations in that the 
inferred magnetic field direction in the plane of the sky 
deviates by an angle of 30° from that expected (parallel 
to the L1544 minor axis). This can be easily explained if 
the cloud is slightly non-axisymmetric (Basu 2000, Galli 
et al. 2001). Thus, while it seems plausible that the mag- 
netic field plays an important role in the evolution of pre- 
protostellar cores just prior to collapse, it is possible that 
the "standard" ambipolar diffusion model requires modifi- 
cation. One notes however that the effects of temperature 
gradients mentioned earlier will have the qualitative effect 
of biasing the 850 /im polarization measurements towards 
the lower density outer parts of cores such as L1544. Thus 
it is clearly of importance to be able to assess how large 
such temperature gradients really are. 

With this in mind, we report here new calculations of 
the expected temperature distribution in pre-protostellar 
cores. We note that a rather similar study has recently 
been carried out by Evans et al. (2001, hereafter ERSM). 
These authors have confined themselves to models with 
spherical symmetry and have compared their model pre- 
dictions to the observations of Shirley et al. (2000). We, 
in contrast, have employed a technique which allows us 
to simulate non spherically symmetric situations such as 
the magnetic field dominated models mentioned earlier. 
We however concur with ERSM in their conclusion that 
the decrease of temperature towards the centers of cores 
strongly affects the observed sub-millimeter emission. We 
also use our models to predict the emission distribution 
expected on the basis of the CB models and conclude that 
a reasonable fit to the observations can be obtained. 

The structure of this paper is as follows. In section 2, 
we describe the model which we have used and the as- 
sumptions which we have made. In section 3, we present 
an analytic treatment aimed at determining the dust tem- 



perature at the center of a spherically symmetric cloud. 
In section 4, we give numerical results for both spherically 
symmetric models and models based on the density distri- 
butions of the CB models. Then, in section 5, we consider 
the intensity distributions which would be expected for 
mm-submm maps of our model cores and compare with 
observed data. In section 6, we summarize our conclusions. 

2. Model 

Our model calculations are aimed at calculating grain 
temperatures under physical conditions similar to those 
thought to pertain in LI 544 and similar pre-protostellar 
cores. We have made a number of simplifying assumptions 
in doing this which we justify in the following discussion. 
One of these is that we neglect scattering as we are in- 
terested in studying cores of more than 10 magnitudes 
of visual extinction which are mainly penetrated by the 
external infrared interstellar radiation field for which the 
Rayleigh limit holds and thus scattering becomes negligi- 
ble. 

Another important assumption which we make is that 
grains are only heated by the incident radiation field and 
that other sources of grain heating are negligible. ERSM 
provide a detailed discussion of the range of applicabil- 
ity of this and conclude that for representative conditions 
with dust grain temperatures above 5 K, heating due to 
the ambient radiation field dominates. We follow them in 
this but differ from them in only considering heating due 
to the (attenuated) external radiation field. That is to say, 
we assume that our model cores are optically thin to their 
own radiation and that we can neglect re-absorption of 
radiation emitted from within our model itself. This is an 
essential simplification which permits us to consider non 
spherically symmetric geometries. However, it restricts us 
to considering structures with hydrogen column densities 
less than a few hundred visual magnitudes of extinction. 

The rationale behind this can be understood if one 
notes that for typical temperatures of order 10 K, one ex- 
pects the bulk of the emission from a pre-protostellar core 
such as L1544 at wavelengths of order 200-300 /xm. The 
available observations (see Andre et al. 2000 for the case 
of L1544) fully confirm this. Then to maintain the optical 
depth below, say, 0.1 at 200 /im requires a column density 
of H2 less than 5 x 10^^ cm^^ corresponding to roughly 
500 magnitudes of visual extinction. Available observa- 
tional data suggest that this condition is satisfied in all 
cases. We note however that for models with a singular- 
ity at the origin such as the singular isothermal sphere of 
Shu (1977), this condition is strictly speaking not fulfilled. 

We determine the grain temperature at a given po- 
sition simply by applying the "classical" equilibrium be- 
tween grain cooling and heating for a spherical grain at 
position r in the dust cloud (see e.g. Spitzer 1978 or 
Boulanger et al. 1998 for a discussion). 



Q.B,[Ti{r)] diy 



QvJu{r) dv. 



(1) 
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Here, the right hand side describes the grain heating due 
to an incident radiation field of average intensity Jiy(r) 
upon grains of absorption efficiency Qi, at frequency v. 
The left hand side gives the cooling rate for a grain of 
temperature Tii{r) and Bi, is the Planck function. 

The only real computational difficulty here is posed 
by the incident radiation field Ju{r) which we suppose to 
be given by the attenuated interstellar radiation field .PJ' . 
Thus, we have : 



Mr) 



ps 



exp[-T^(r, 6*, 0)] dVt 



(2) 



Here, the integral over solid angle describes the average 
attenuation due to an optical depth Ti,{r, 6, (p) in any given 
direction. 

In order to compute grain temperature for a model 
cloud of relatively arbitrary geometry, one thus must find 
an efficient method of performing the angle integration 
of equation 2. Once the grain temperature is determined 
for all positions, a simple integral allows one to derive the 
intensity distribution expected at a given wavelength. One 
merely needs to choose the incident interstellar radiation 
field J]^ and the grain absorption efficiency Q^. 

There are a number of choices which one can make con- 
cerning the grain opacities. Mathis (1990) has tabulated 
values expected for "standard" interstellar grains and we 
have used these for comparison purposes. However, within 
dense dust clouds, one expects grains to acquire ice man- 
tles and this will substantially change their optical proper- 
ties. In addition, one expects at the densities of interest to 
us that grain coagulation may change the size distribution. 
We have therefore used for many purposes the opacities 
computed by Ossenkopf & Henning (1994, hereafter OH) 
for grains which have coagulated for a time of 10^ years. 
As a standard case, we have used their case for thin ice 
mantles and density 10® cni^'^ but we note that in many 
situations, it is more realistic for pre-protostellar cores 
to consider their "thick ice mantle case". We have found 
however that for our purposes, the difference between the 
two is not great as is the difference between their solu- 
tions for different densities. This is consistent with the 
fact that Kramer et al. (1997) find that the observed ratio 
of infrared extinction to millimeter emissivity in IC5146 
(for visual extinctions in the range 20-30 magnitudes) is 
compatible with the OH results. 

In similar fashion, we have used as a comparison stan- 
dard the solar neighbourhood radiation field given by 
MMP. Black (1994) has updated this work using results 
from the FIRAS experiment on COBE (Wright et al. 
1991). This suggests that the MMP field is an underes- 
timate in the mid and far IR . However, in any case, 
the values used are an "educated" guess concerning the 
field actually incident upon cores such as LI 544 and one 
can expect differences from one region to another. For ex- 
ample, it is clear that small particle emission in the mid 
infrared can be an important component of the incident 
radiation field and one expects this to depend on the far 
ultra-violet radiation from neighbouring 0-B stars. This 



can be expected to vary depending on how many mas- 
sive stars have formed in the vicinity of a given core. One 
notes from the work of Bacmann et al. (2000, their table 
1) that the intensity in the LW2 ISOCAM filter (S-S.S/xm) 
is 0.6 MJy/sr towards L1544 in Taurus (a factor of order 3 
larger than the average radiation field at this wavelength 
according to Black (1994)) but 20.6 MJy/sr towards the 
OphD core (after correcting for zodiacal light). This seems 
likely to be typical also of the small particle radiation at 
wavelengths up to 50/xm. 

Finally, we need to define the density distribution of 
our models. For spherically symmetric cases, we have 
used the following density distributions: (z) a homoge- 
neous sphere, (ii) a Bonnor-Ebert sphere (Bonnor 1956, 
Ebert 1955), and [iii) a singular isothermal sphere 
(Chandrasekhar 1939, Shu 1977). Case {i) was chosen 
to test the accuracy of the numerical scheme, since in 
this simple case exact solutions for the mean intensity 
and analytical approximations for the dust temperature 
are easily derived. Cases {ii) and {in) represent equilib- 
ria of self-gravitating non-magnetized isothermal clouds 
and have been extensively adopted to model the density 
profiles of pre-protostellar cores. Notice that in cases {ii) 
and {iii) the structure of the cloud is computed assum- 
ing that the gas is isothermal. However, the dust clearly is 
not. Moreover in the central parts of the core at densities 
above 10^ cm^^, the gas and dust temperatures are likely 
to be coupled (see e.g. Goldsmith 2001) and hence there 
is strictly speaking an inconsistency which we will ignore 
in the following discussion. It is however a point which 
requires further study. 

In all cases (including the non spherically symmetric 
models discussed below), the model "core" was truncated 
at an outer radius i?out, and embedded in an envelope of 
radius i?onv and uniform density ponv, equal to the density 
of the cloud at radius i?out- The rationale of this is that 
it simulates the "outer skin" of the molecular cloud where 
a PDR (Photon dominated region) shields the molecular 
gas against UV and is consistent with the idea that the 
internal molecular cloud pressure should not fall below the 
average interstellar pressure. Thus, the value of penv was 
chosen to simulate the effect of the internal CMC pressure 
Pg confining the core. One expects Pq to be a factor of a 
few larger than the mean ISM pressure of roughly 2 x 10^ 
K cm^^(sec section 2 of McKee 1999). If we take Pg equal 
to this minimum value and assume a typical line width of 1 
kms~ we find an envelope density corresponding to riH2 
— 1000 cm^'^. This outer PDR layer should provide shield- 
ing equivalent to 1 magnitude of visual extinction and we 
define _Ronv such that riHa (-Renv — -Rout) = 1 X 10^^ cm^^. 
We note that while these assumptions are somewhat ad 
hoc, they have little influence on our results which con- 
cern the high density core where the extinction is much 
higher than 1 magnitude. 

Our non spherically symmetric models are based on 
the results of CB who have studied the development of in- 
fall driven by ambipolar diffusion assuming axial symme- 
try. They present results for a "disk-like" model and give 
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the mid-plane density as a function of radius at times tj 
{j = 0,1,2,3...) when the central density nc{tj) = 10-' tIco 
with nco = 4.4 x 10^ cm"^. The times tj are 0, 2.27, 
2.60, 2.66, 2.68 ... Myr from which one sees that the 
last stages of the evolution are extremely rapid. In this 
study, we consider the density distributions at times ^2, 
^3, and ti as being of most interest from the standpoint 
of comparison with data for cores similar to L1544. We 
also need to make an assumption about the density dis- 
tribution in the direction perpendicular to the disk mid- 
plane. We have assumed that hydrostatic equilibrium is 
valid in the z direction (perpendicular to the mid-plane) 
and hence that the density can be expressed as n(r^ z) = 
nc{tj)FcB{'i')4>{z) with Fcsir) taken from the study of 
CB and 0(z) = l/cosh^(z/_ff) where the scale height is 
H = [kT/{2T:Gp)]^/^ and p is the density. For this pur- 
pose, the disk is considered to be isothermal at a temper- 
ature of 12 K. Again, this is certainly not self-consistent 
for the reasons mentioned earlier but this seems a reason- 
able approximation to the density distribution of a CB 
model. 

3. Analytical results for spherically symmetric 
clouds 

The results for spherically symmetric models are impor- 
tant to us both because they can be treated to some extent 
with analytic techniques and because direct comparisons 
are possible with ERSM. As a check on our computational 
accuracy, we first show our results for the mean intensity 
in a homogeneous spherical cloud compared with the an- 
alytic solution. 

3.1. Mean intensity in a homogeneous sphere 

The solution of the full equation of radiative transfer in 
a homogeneous sphere can be obtained analytically with 
the spherical harmonics method (Flannery, Roberge, & 
Rybicki 1980). In Appendix A we show that for a homo- 
geneous cloud of radius R and centre-to-edge optical depth 
Ti/c, the calculation of J,/{r) from eq. (0) is straightforward 
and can be expressed in terms of exponential- integral func- 
tions of the argument x± — t^c(1 ± x), where x = r/ R: 



2{x+ — X-) 
[x+E2{x^) - x_E2{x+) + x_Eo{x_) - x+Eo{x+)] . (3) 

Approximate expressions for the mean intensity can be 
obtained in the limit of very large or very small t^c- In the 
optically thin hmit eq. (0) becomes 



(4) 



j,{x)^.rju 



1 l-x\ 


n + x\] 


- H In 





2 4x 


^l-xj\ 



whereas an approximate formula (diverging for x —>■ 1) for 
the optically thick case is 



0.96 



0.94 



0.92 - 



0.9 




O.S 0.4 0.6 0.8 1 
r/R 



0.2 0.4 0.6 0.8 1 
r/R 



J^{x) ~ J, 



^ ' l/C "^ \ 



(5) 



Fig. 1. Mean intensity Ji, for a homogeneous spherical 
cloud as a function of radius x = r/R computed by nu- 
merical integration of eq. (0) {solid lines) compared with 
the analytic formula eq. (0) {triangles). The approximated 
expressions (0) and (ph for the optically thin and thick 
cases are shown by dashed lines. Results are shown for 
two model clouds with Ti,c. — 0.1 and 10. 



In Figure nl we compare the exact and approximate 
analytical results with that obtained numerically integrat- 
ing eq. (ph. We see that the numerical code reproduces the 
analytic result very well. Thus, in simple geometries, the 
procedure used for the angle integration is adequate. 

3.2. The dust temperature at the centre of a spherical 
cloud 

There are great advantages in developing analytic formu- 
lae which allow an approximate estimate of the tempera- 
ture in dense pre-protostellar cores. This permits for ex- 
ample a discussion of scaling relationships for different ex- 
ternal radiation fields and opacities. We have accordingly 
attempted to derive an analytical estimate of the temper- 
ature at the center of a spherically symmetric cloud for 
comparison with our numerical results. 

In order to develop an analytic approximation, we 
must first make some crude approximations concerning 
the dust opacity and external radiation field in cores such 
as L1544. Following Black (1994) and MMP, the inter- 
stellar radiation field at wavelengths above 0.1 /Lim can 
be represented by the sum of four contributions: (1) an 
optical-NIR component peaking at Ap = 1 /im, due to the 
emission of disk dwarf and giant stars; (2) the diffuse FIR 
emission from dust grains, peaking at Ap = 100 //m; (3) 
mid-IR radiation from small non-thermally heated grains 
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in the range 5-100 /im and (4) the cosmic background 
radiation, peaking at Ap = 1 mm. Components 1,2, and 
4 can be represented by a single (or a sum of) modified 
black-bodies at temperatures Ti of the kind 



J'f = 



2hv^ fXpV 

T 



w, 



ohu/kTi _ I ■ 



(6) 



where Ap is the peak wavelength and Wi are dilution fac- 
tors. Th e val ues of the parameters Ap, Wi and Ti are listed 
in Table BJ. The third (MIR) component is found to be 



variable near dense cores such as L1544 and can be ap- 
proximated as a power law. In the following discussion 
however, we at first neglect it and then make an estimate 
of its importance. 

As in the numerical computations, we have based our 
opacity law on the work of OH (thin ice, density 10^ 
cm^^). For the purposes of our analytic work, a very crude 
piecewise power law fit with breaks at 10 fim and 400 /im 
has been used (see Table B.2). We then assume that in 
a sufficiently large range of wavelengths around the peak 
of each component, the dust absorption efficiency Q^ can 
be approximated by a power-law, Q^, oc (Ap/A)". The val- 
ues of the parameters Ap and a used by us are listed in 
Table 



B.2 



Let the radiation absorbed (indicated by the subscript 
"a") by dust be concentrated around a frequency i>a, and 
the radiation emitted (indicated by the subscript "e") by 
dust be characterized by a frequency Ve- We take Ve and 
Qfe equal to the values for absorption in the FIR (see 
Table B.2), and we consider the heating of the dust by 



each component of the radiation field separately. 

With these assumptions, the mean intensity at the 
cloud's center is 



J^(0) = fj exp {-T^c) = J'J' cxp 



(7) 



where Ta is the centre-to-edge optical depth of the cloud at 
frequency Va- Then, defining Pd = kTd/hve-, Pi = kTi/hva, 
eq. (m) becomes 



4 / \ p 



^=f(-j (^ji:»'A> (8) 



where Qa and Qe are the values of Qy at frequencies Va 
and i^e, respectively, and 



/■oo j-ae+3 



-^7^^-^ dt = r(ae + 4)C(«e + ^)PT^\ (9) 



B, 



f 



at/f}. 



1 



exp (-Tat"") dt, 



(10) 



with q = tta + P + 3. In the above, F and C, are the 
gamma and Riemann zeta functions (see Abramovitz & 
Stegun 1965). The integral in the equation for Bi can be 
easily evaluated in the two limiting situations r^ ^ 1 and 



Tq ^ 1 . For Ta <C 1 the exponential in Bi can be expanded 
as a Taylor series, giving, to first order in Ta, 



B, ^ r{q + l)C{q + l)Pf+' 
-TaT{q + aa + l)C{q + a 



m 



q+aa+l 



(11) 



Thus, the dust temperature at the center is given, to 
lowest order, by 



kTd 
hv. 



-±^c\y mp: 



9-1-1 






where 



C = 



r(g + i)C(g + i) Qa 

F(ae -f 4)C(ae -I- 4) Qe 



(12) 



(13) 



We note here that for the case of interest in which the 
dust is heated by (optically thin) FIR radiation, Qa = Qe, 
i^a = Vp = ^e and one finds that C ~ 1.68 is a constant 
slightly dependent on the power law index of the dust 
opacity but otherwise independent of dust characteristics. 

For Ta ^ 1 the integrand in Bi can be expanded as a 
Taylor series for small i, obtaining the approximate result 



B, 



^P^_ 

aa 

1 

2Q:a 



-F 



--q/aa 



-(q+l)/aa 



+ 



(14) 



Thus, the dust temperature at the center is given, to low- 
est order, by 



^ _ n^-g/[("e+4)aa] 

7 — a 

hve 
where 



D 



r(— ) 



; w,p^ 



Qa 



aaT{ae + 4)C(ae + 4) 



(15) 



.(16) 



Notice that the dependence of the dust temperature on 
the intensity of the external radiation is rather weak, Td oc 
l^i/(Qe+4) Again, this expression simplifies considerably 
when the dust is heated by (optically thick) FIR radiation, 
and D ~ 0.59 in this case. 

The above discussion has neglected the contribution 
of the MIR radiation. We approximate the external MIR 
field with a power law spectrum in the wavelength range 
10-100 /im and with a lower frequency cut-off at 100/xm, 



2hvl 



(17) 



with the values of W , Ap and p given in Table B.l , Notice 
that p < in this case and thus for the "MIR" radiation, 
there is approximate cancellation in the integral on the 
right hand side of equation |l| between the frequency de- 
pendence of Qi, and that of the external MIR field. The 
consequence is that the main contribution to the integral 
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comes from frequencies where the optical depth is of or- 
der unity (i.e. at a wavelength of roughly 20/im for 74^=30 
mag.). 

Substituting the above expression for the MIR field in 
eq. (|l]), we obtain 

4 






r(s,TQ) 



aa 



(18) 



where s = (aa + p + l)/aa = 0.5 for our parameters, 
and r(s, Ta) is the incomplete gamma function. Here Tq is 
the cloud's optical depth at 100 /xm, which becomes unity 
at Av — 414. Below this value, the incomplete gamma 
function can be expanded for small r^, obtaining 

s+l 

+ .... (19) 



r(s,r,)^r(s) + ^-^ 

s s + l 

In eq. ([l^), ([l5|), ( [l8|) and (|l9|), it is convenient to ex- 
press the optical depth Ta in terms of the optical extinction 

Av = IMQtv = 1.086 (^\ Ta, (20) 

where the ratio in parentheses can be comput ed fo r each 
components with the parameters listed in Table [R^ for the 
Ossenkopf & Henning (1994) opacity. We see that a cloud 
is optically thin to V-NIR radiation for Av ^ 3, to FIR 
radiation for Av <C 700, and to the cosmic background 
radiation for Ay < 6 x 10"^. Thus, with Black (1994) 
standard interstellar field, the central temperature in a 
"diffuse" cloud core of typical centre-to-edge extinction 
Av — 1-2 is determined by optically thin absorption of 
V-NIR and MIR radiation (the effects of an external UV 
field are not considered here) whereas for 10 < Av ^ 400 
is determined by optically thick absorption of V-NIR ra- 
diation and optically thin absorption of MIR and FIR ra- 
diation. 

Inserting the numerical values given in Appendix B 
(and keeping terms up to second order in Ta), the formulae 
above give 



Trf(V-NIR, thin) ~ 16.5 - 3.6A; 
Td (V-NIR, thick) ~ 43 A 



-0.56 



77 Ay^-^^ 



(21) 
(22) 
(23) 



rd(FIR, thin) ~ 6.2 - 0.003lAy + . . . , 

Td(MIR, thin) ~ 7.9Ay°°^^ 

X (1.8 - 0.098^^^ + 7.9 X W'^Al}^ + . . .)i/5^ (24) 

where Td is in K. 

Combining eq. (El])-(p3), one can obtain simple ana- 
lytical expressions for the dust temperature (in K) at the 
cloud's centre, 

Td ~ [Td(V-NIR, thin)5-6 + T^MIR, thin)^-6]i/5.6^ (25) 
for Av — 1-2, and 

Td ~ [Td(V-NIR, thick)^-^ -I- rrf(MIR, thin)^-^ 

-|-rd(FIR, thin)5-6]i/5-^ (26) 




50 100 

visual extinction Ay 



150 



Fig. 2. Dust temperature Td at the centre of spherical 
clouds of visual extinction Av- The contribution of V- 
NIR radiation (dotted line), MIR radiation [long- dashed 
line), FIR radiation (short-dashed line), and the total 
(thick solid line), obtained with the analytical formulae 
of Sect. 3.2 are compared with the numerical results (tri- 
angles). The inset shows a blow-up of the figure for ex- 
tinctions less than 11 magnitudes. 



for 10 < Av ^ 400. Notice that increasing the intensity of 
the radiation field by a factor g results in an increase of 
the dust temperature by a factor g^^^-^. 

We have checked these approximate formulae against 
the numerical computations and the results are shown in 
Figure ^ Here we compare in the first place the numerical 
model predictions for the temperature at the center of a 
Bonnor-Ebert sphere with that derived from eq. (|2^) and 
(P6|). We see that the analytic formula is in agreement 
with the more accurate numerical results within ~ 1 K . 

It is also of interest to compare the expectations from 
the analytic formula with that calculated as a function of 
radius in a given Bonnor-Ebert sphere. Here, we merely 
use the "intuitive approximation" that the appropriate 
optical depth to substitute in eq. (|2^) is the minimum op- 
tical depth along a ray to the exterior of the core. We see 
in Fig. H that this approximation is also quite reasonable 
suggesting that, in a number of situations, one can obtain 
an estimate of the dust temperature as a function of posi- 
tion using eq. (p6|). Obviously however, this must be used 
with caution. Nevertheless, it should as a rule be more 
accurate than the isothermal assumption. 
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the BES distribution "flattens out" . The computed tem- 
perature distributions for the two cases are close to being 
identical even in the central regions. There is in both cases 
a substantial fall off from values above 15 K at the edge (in 
the constant density envelope mentioned above) to values 
as low as 5-6 K at a radius of 10^'^ pc or 200 AU. 

These results do not change greatly for different choices 
of grain opacity and radiation field. However, a larger ra- 
diation field, as estimated by Black (1994) in the range 
10-50 /im, causes a slight increase in the central dust tem- 
perature. In general, we conclude that in structures such 
as these (with visual extinctions of order 100 mag.), the 
central temperature can be expected to drop by roughly a 
factor of 2-3 relative to that measured on scales of ~ 0.1 
pc. We have compared our results with those of ERSM 
for similar parameters and find that they obtain slightly 
higher temperatures (7.5 K as compared to 7.2 K for the 
same conditions) which may be due to our neglect of heat- 
ing due to the radiation of the core itself. However, the 
general behavior seems very similar for both codes al- 
though, as noted above, our neglect of heating by core 
radiation is not valid for the SIS distribution at small radii. 
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Fig. 3. Comparison of the radial dependence of dust tem- 
perature computed for model Bonnor-Ebert spheres with 
central extinctions 30, 90, and 150 magnitudes {dotted 
curves) with results from analytic approximations (thick 
full lines) in eq. ( |2^ ) and (p5|). Here, the extinction Ay 
in the analytic formulae has been interpreted as the min- 
imum extinction along a ray emerging from the core. 



4. Numerical computations of the dust 
temperature distribution 

In this section, we present our results obtained numerically 
integrating eq. @) and (g). 



4.1. SIS and Bonnor-Ebert spheres 

We first consider inhomogeneous spherically symmetric 
density distributions. Two obvious cases of interest are the 
singular isothermal sphere (SIS) solution and the Bonnor- 
Ebert sphere (BES). In Fig. ^, we show the results we have 
obtained for these density distributions. Both models are 
obtained from the equation of hydrostatic equilibrium as- 
suming a gas temperature of 12 K. For the BES, we have 
chosen a model with a central gas density riHa = 4.4 x 10^ 
cm~^ similar to some observed pre-protostellar cores. One 
notes first that the two density distributions only differ 
greatly at small radii (less than 10^ AU or 0.03 pc) where 



4.2. Temperature distribution in the CB model 

We now consider the temperature distribution to be ex- 
pected within the model proposed by CB for L1544. As 
discussed in section 2, this is a disk-like axisymmetric 
model. In this discussion, we consider the radial distri- 
bution of dust temperature in the mid-plane of the disk 
as well as the temperature distribution along the axis of 
symmetry perpendicular to the disk. The calculations are 
presented for the density distributions corresponding to 
three evolutionary times ^2 = 2.60 Myr, t^ = 2.66 Myr, 
and ^4 — 2.68 Myr during which time the central density 
evolves from a molecular hydrogen density of 4.4 x 10^ 
cm~'^ to 4.4 X 10^ cm~'^(see section 2). 

In Fig. we present results for these models. One notes 
the fact that the CB models are relatively "thick" disks in 
that their extent along the z-axis is an appreciable frac- 
tion of their extent in r. The dust temperature computed 
by our models behaves in similar fashion to that obtained 
in our spherically symmetric models. The main effect of 
the disk structure is simply that one finds dust at temper- 
atures of 10 K or higher much closer to the center of the 
core than in the spherically symmetric case. 

Another view of the temperature variation is presented 
in Fig. H which shows a cut perpendicular to the disk mid- 
plane. One sees that close to the midplane, the attenua- 
tion of the external radiation field is large and the tem- 
perature remains relatively low while in the perpendicular 
direction, it rises more rapidly. The temperature gradient 
is larger along the minor axis if one examines this model 
edge on. In practise, one must take into account effects 
due to finite inclination as well as the contributions from 
foreground and background layers. 
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Fig. 4. Dust temperature as function of radius computed 
by the numerical model for the SIS [jull lines) and Bonnor- 
Ebert sphere (central density 4.4 x 10^ cm^"^, dashes) dis- 
tributions. The top panel shows the density distribution 
and the bottom two panels show the temperature distri- 
bution we compute for two different choices of external 
radiation field and grain optical parameters. The center 
panel shows results for the MMP radiation field together 
with the grain opacities recommended by Mathis (1990). 
The bottom panel shows results obtained using the opac- 
ities of OH (their col. 5) and the radiation field of Black 
(1994). 



5. Comparison of Models with Data 

In this section, we present calculations of the mm-submm 
emission predicted by the models discussed above. We 
compute the expected intensity distributions and compare 
with observations of L1544. We consider first the spheri- 
cally symmetric cases. 

5.1. Intensity Distribution for spiierically symmetric 
models 

Using the dust temperature distribution discussed earlier, 
we have computed the expected intensity distributions at 
three wavelengths for which observational data are avail- 



Fig. 5. Dust temperature variation with position in the 
models of CB. We show in the top panels the density dis- 
tribution implied by their model in the midplane of the 
disk (top left) as a function of radius r and perpendicular 
to the disk as a function of height above the midplane z 
(top right). Results are shown for three different times in 
the evolution: at 2.60 Myr (^2) when the central II2 den- 
sity is 4.4 X 10^ cm"-^, at 2.66 Myr (^3) when the central 
density is 4.4 x 10^ cm~^, and at 2.68 Myr when the cen- 
tral density is 4.4 x 10^ cm~'^. The bottom panels show 
the corresponding temperature distributions computed by 
our code as a function of r and z. Dashes are for time ^2, 
the full line for ^3, and dots for t^. 



able in the literature. These results have been obtained 
simply integrating the transport equation (with our calcu- 
lated temperature distribution) across our model core. We 
have convolved the resultant intensity distribution with a 
"gaussian beam" of 13, 15 and 8" (at 140 pc) for 1.3 mm, 
850, and 450 iiui, respectively. We present in Fig. re- 
sults using the Black (1994) radiation field and compare 
the expected intensity profile at wavelengths mentioned 
above from the BES and SIS models. 

One sees that at all wavelengths but especially at 450 
/im, the profile is much flatter using our calculated tem- 
perature distribution than in the isothermal case. This 
just reflects the fact that the cold core nucleus predicted 
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Fig. 6. Dust density and temperature variation in the 
model of CB at time ts (2.66 Myr with central density 
4.4 X 10^ cm~^)seen in the plane perpendicular to the disk 
midplane. The density variation shown in color is that of 
a flared disk. The temperature variation is given by the 
contours. 



by our calculations does not contribute substantially at 
the shorter wavelengths. Even at 1300 /zm, the contrast 
between the intensity at offset 10~^ pc (200 AU) and 0.1 
pc is 26 in the isothermal case and 9 using our computed 
temperature gradient. Thus we conclude (in agreement 
with ERSM) that conclusions about the density structure 
based upon the intensity distribution of the dust contin- 
uum emission must take temperature gradients into ac- 
count. 

Another clear result from Fig. M is that the predicted 
intensity distributions for the BES and SIS models are 
essentially identical. The low temperature of the central 
nucleus causes the observed intensity profile to be insensi- 
tive to the density distribution inside 0.01 pc (2000 AU). 
This means that one should be very cautious about claims 
for a flat central density distribution based upon maps of 
millimeter dust emission. On the other hand, derivations 
of the density profile based upon either NIR reddening 
(Alves et al. 2001) or mid-IR extinction (Bacmann et al. 



Fig. 7. Intensity as a function of projected radius for the 
SIS and BES models discussed in section 2 and shown 
in Fig. 0. The top panel shows results at a wavelength 
of 1.3 mm, the center panel at a wavelength of 850 /im 
and the bottom panel at a wavelength of 450 ^m. We 
show results for the Bonnor-Ebert sphere {dashed) and 
SIS models (dotted line) with the Black radiation field and 
OH opacities. For comparison, we also show the expected 
distributions for a SIS model at a dust temperature of 
12 K {soM). 

2000) are insensitive to dust temperature and hence are 
to be preferred for such purposes. 

5.2. Intensity Distribution for the CB model 

The spherically symmetric models discussed above clearly 
do not apply to cores such as L1544 which has an ex- 
tremely elongated intensity distribution when mapped at 
millimeter wavelengths (as well as an elongated appear- 
ance when observed in absorption at 7 /xm , Bacmann et 
al. 2000). We therefore now consider the expected inten- 
sity distributions for the models of CB shown in Fig. ^ 
One should note that these were originally constructed to 
fit the observed millimeter continuum intensity distribu- 
tion for LI 544 assuming isothermal dust. We now con- 
sider the consequences of the more realistic temperature 
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Fig. 8. Intensity distribution {full lines) for wavelengths 
of 1.3 mm (top), 850 /xm (center), and 450 fim (bottom) 
predicted by the CB model along the major axis (left hand 
panels) of L1544 and along the minor axis (right hand 
panels) assuming the density distribution at time is (2.66 
Myr) with central density 4.4 x 10® cm^'^ . In the panels 
for 850 /xm we also show profiles computed for the CB 
model at time ^4 {dash-dot, 2.68 Myr) and ^2 {dots, 2.60 
Myr) (see caption for Fig. 0). An inclination of 16° to the 
plane of the sky has been assumed. For comparison, we 
show also observed cuts along the major and minor axes of 
L1544 taken from Ward-Thompson et al. (1999) (1.3 mm 
data) and Shirley et al. (2000) (850/im and 450 /im data). 
The model intensity distribution has been convolved with 
gaussian beams with an HPBW of 13, 15 and 8 arc seconds 
for 1.3 mm, 850 /im and 450/im, respectively, i.e. the same 
observational beam as the data. 



distribution discussed above and compare with the avail- 
able observations. We have again here used dust opacities 
from OH (thin ice, col. 5) and the external radiation field 
of Black (1994). 

In Fig. y, we show the resulting intensity distributions 
for the epoch is of CB (central density 4.4 x 10® cm~^). 
We have assumed an inclination of 16° with respect to 
the plane of the sky as suggested by CB and compare 
with cuts along the major and minor axes from the work 



;=37° and Th=12 °K 




wavelength [/J-nn] 



Fig. 9. Observed spectral energy distribution for L1544 
compared with model predictions {solid lines) for the 
flux in a 40" aperture as a function of wavelength. The 
model results are for the flux (Jy) in an aperture of 40" 
HPBW towards the peak of the L1544 dust emission. 
Observational data (error bars) are from Shirley et al. 
(2000) and Ward-Thompson et al. (1999). The model data 
are for 3 different inclinations of a CB disk (time is) to 
the plane of the sky. The dashed line shows for compari- 
son the expected SED for an isothermal CB disk inclined 
at 37°. 



of Ward-Thompson et al. (1999) and Shirley et al. (2000). 
We conclude that despite departures from isothermality, 
the agreement between model and observed intensity dis- 
tribution is reasonable at all three frequencies for time t^ ■ 
It certainly gives a better fit than the profile for time ^2 
although it would be difficult to distinguish between t^ 
(2.66 Myr) and t^ (2.68 Myr) on this basis. Thus we con- 
firm the result of CB that this stage in the development 
of their model fits the observed intensity distribution. We 
note nonetheless that in all cuts along the major axis, 
the fall off in observed intensity at positive offsets is more 
rapid than in the models. 

Another characteristic that a model should fit is the 
observed spectral energy distribution or SED. A com- 
plication here is that observations typically resolve pre- 
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protostellar cores and that the SED thus depends sensi- 
tively on the assumed aperture. In Fig. 0, we present the 
expected spectrum for L1544 compared with observations 
from Shirley et al. (2000) for an assumed 40" aperture (a 
compromise value). One sees that edge-on disks of this 
type are considerably stronger and easier to observe than 
face-on disks. This is because (even with a 40 "beam at a 
distance of 140 pc) we resolve the source and as at these 
wavelengths, the model is optically thin, the flux is greater 
for the edge-on case where the column density is larger. 
We show for comparison the SED for the same model den- 
sity distribution assumed to be isothermal and with an in- 
clination of 37° . The isothermal model is slightly "hotter" 
than the model with computed temperature distribution 
and has a somewhat steeper long wavelength spectrum. 

We get a reasonable agreement with the observed SED 
for an inclination of 37° though this result reflects the as- 
sumed dust opacity and should be treated with caution. 
There is still considerable uncertainty concerning the na- 
ture of the ice mantles on the dust grains at the heart of 
pre-protostellar cores like L1544. We conclude neverthe- 
less that the non-isothermal CB model at time t^ does 
indeed give a reasonable fit to the observed flux distri- 
bution. The effects of departures from isothermality on 
the SED are probably small in comparison to the observa- 
tional difficulties (e.g. due to pointing errors) of obtaining 
matched apertures at different frequencies towards an ex- 
tended source. This is both because the fluxes from a 40" 
region are not very sensitive to the central density distri- 
bution (i.e. the central 2000 AU) and because the low tem- 
perature causes a high density nucleus (if present) to have 
much less weight. This is in particular true at 850 /im and 
higher frequencies and from this point of view, 1.3 mm ob- 
servations of high angular resolution are needed to probe 
the nucleus of such cores. 

As seen above, if we interpret cores like L1544 in terms 
of the CB model, the inclination of the CB disk to the line 
of sight is a critical parameter. It is also important for the 
interpretation of Zeeman effect observations (Crutcher & 
Troland 2001). This suggests that, for a model like that 
of CB, it is useful to estimate the inclination based upon 
the ratio of major and minor axes in a tracer such as 
the millimeter continuum. This must be done using so far 
as possible the true intensity distribution of the model 
core. With this in mind, we have attempted to calibrate 
the measurement using CB's model at time t^ and show 
the results in Fig. ^ We compare the results from the 
numerical computations with the formula given by CB for 
a disk of finite thickness as well as the simple result for 
a thin disk. Our computed axis ratio behaves rather like 
a thin disk at moderate inclinations but is close to the 
expected result in the edge-on case. We obtain from the 
axis ratios a best fit for an inclination of 37° consistent 
with the SED of Fig. ^ shown above. However, our data 
would certainly accomodate an inclination of 25° and (as 
Fig. H demonstrates) we cannot exclude 16°. 
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Fig. 10. Ratio of half-power sizes along the major and 
minor axes computed for CB model t^ as a function of in- 
clination of the mid-plane relative to the line of sight (full 
line). The measured ratio for L1544 is shown for compar- 
ison (dashed horizontal line) with upper and lower limits. 
The dotted lines show the axis ratios expected for a ge- 
ometrically thin disk (sini) and for a geometrically thick 
disk using eq. (1) from CB (upper dotted curve) with ratio 
of disk thickness to radius 0.27. 



6. Discussion and Conclusions 

This study has shown that departures from isothermal- 
ity in pre-protostellar cores, while not a dominant effect, 
should be taken into account when considering core struc- 
ture at times just prior to the onset of collapse. This is of 
importance when attempting to deduce the density struc- 
ture of such regions based on maps of dust emission. It 
also may be important when considering the stability of 
such cores. Above densities of 10^ cm~^ , one expects gas 
and dust temperatures to be coupled (see Goldsmith 2001, 
Kriigel & Walmsley 1984). Thus at higher densities in the 
center of cores such as L1544, one can expect the gas tem- 
perature to follow the dust temperature. A decrease in 
dust temperature in the nucleus of such cores requires a 
steeper density gradient than in an isothermal model and 
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may have some consequences for the behavior of the mass 
accretion rate with time. 

We have confirmed the ERSM conclusion that it is 
essentiahy impossible to distinguish a highly peaked den- 
sity distribution (such as an SIS) from one with a flat- 
tened center (such as BE spheres) based on the millimeter 
maps. This has of course no consequences for methods of 
determining the dust column density distribution based 
upon absorption or extinction measurements (e.g. Alves 
et al. 2001, Bacmann et al. 2000). These other approaches 
however have their own limitations. For example, using 
NIR colors one is in practise limited to background gi- 
ants seen in both H and K which puts an upper limit on 
the possible extinctions which arc probed. The MIR tech- 
nique on the other hand has problems both because one 
needs to distinguish between foreground and background 
radiation fields and because one is forced to make assump- 
tions about the uniformity of the background field. Both 
these effects make it difficult to provide incontrovertible 
evidence of extinctions of order or larger than 100 mag- 
nitudes if present. We believe therefore that one should 
still be very cautious about ruling out peaked density dis- 
tributions (such as the SIS) on the basis of the available 
data. 

Using for the first time a reliable dust temperature 
distribution, we have obtained reasonable agreement be- 
tween our models based upon the CB density distribution 
at time t^ corresponding to a central density of 4 x 10^ 
cm~^ and the available millimeter and sub-millimeter data 
for L1544. This is a strong argument in favor of the va- 
lidity of the model of CB as far as L1544 is concerned. 
We favor however higher inclinations of L1544 to the line 
of sight than given by CB and can fit both the SED and 
the aspect ratio for an inclination of 37°. There is on the 
other hand a considerable margin of error in this estimate 
and we do not entirely exclude the value of 16° preferred 
byCB. 

We conclude in any case that it would be useful to 
fit other cores with similar models. There may be some 
possibility in this way of relating the observed structure 
to "age" though it remains to be seen how effective this 
is in practise. The high central column density is the best 
indicator of an age close to the pivotal situation where dy- 
namical collapse can commence. Unfortunately, the lower 
temperature of the high density nucleus caused by the in- 
creased extinction makes this high central column difficult 
to discern. 

It should be realised also that other models are likely 
to be capable of fitting observational constraints. Recent 
relevant theoretical studies include those of Safier, McKee, 
and Stabler (1997), GaUi et al. (2001) , and Indebetouw 
and Zweibel(2000). We note that GaUi et al. (2001) ob- 
tain an excellent fit to the observed intensity distribution 
at 1.3 mm of L1544 with their model of a slightly in- 
clined "singular isothermal disk" or SID. This is a non- 
axisymmetric structure with field lines running approx- 
imately in the line of sight. It would be useful also to 
compute the temperature structure of such a model in 



analogous fashion to that used in this paper for the CB 
models. 

We have not considered in detail in this study the effect 
of external radiation fields other than the "standard" fields 
of MMP and Black (1994). However, it is hkely that the 
actual fields incident upon dense cores vary from core to 
core and still more between molecular cloud and molecu- 
lar cloud. Such differences plausibly occur due to the pres- 
ence in the vicinity of nearby OB stars. The Bacmann et 
al. (2000) study (see their table 1 and subtract the zodia- 
cal contribution) suggests that there may be as much as an 
order of magnitude difference in the ISO LW2 (5-8.5 /im) 
background towards different cores. It in fact also suggests 
that we underestimate by a factor of order 3 the field (at 
7 /im) incident on L1544 using the Black (1994) average 
field. We have examined IRAS maps of LI 544 and aver- 
aged intensities in an annulus surrounding the core with 
inner radius 100" and outer radius 300" centred on the mm 
dust emission peak. We conclude that the field adjacent 
to the core is a factor 2.5-4 greater than the Black (1994) 
field at 12 and 25 /im suggesting that we underestimate 
the MIR contribution to the incident field (section 3). This 
may cause our temperature estimate to be low (by a fac- 
tor of roughly 1.2) at extinctions of 20-30 magnitudes (see 
Fig. i. 

We may also be under-estimating the FIR field. L1544 
is quite clearly present on IRAS lOO/im maps though not 
present at shorther wavelengths. This suggests that the 
incident UV-optical radiation field (presumably also re- 
sponsible for exciting the MIR) is converted into far IR 
radiation in a "FDR" layer surrounding the cold core. This 
presumably should be added to the average incident field 
which we have used in our computations. We thus may 
also be underestimating the FIR heating. 

We expect temperatures in these cores to vary as x^^^'^ 
where x is the scaling factor for the external radiation 
field (see section 3) and thus we expect differences as large 
as 1.5 between different objects. We would expect for ex- 
ample the nucleus of a core like OphD on this basis to be 
roughly 1.5 times hotter than L1544 using the mid IR in- 
tensities of Bacmann et al. (2000). This will have some ef- 
fect on parameters such as the accretion rate but probably 
not an important one since the region where gas and dust 
temperatures are coupled is very limited. On the other 
hand, phenomena such as depletion of molecular species 
on grain surfaces are extremely sensitive to grain temper- 
ature and these are likely to have different characteristics 
depending on the precise grain temperature. This may ex- 
plain why some of the chemical characteristics of cores in 
Ophiuchus and Taurus differ (long chain carbon species 
for example are much less prevalent in Ophiuchus) and 
may have consequences for the ionization degree at high 
densities. 

One of the most useful products of this study is the 
analytic estimate which we have made of the temperature 
at the center of a spherically symmetric core. We have 
shown that the formula which we have derived gives in 
many cases a reasonable fit to the radial dependence of 
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dust temperature and we believe it can be used as a useful 
first approximation to the gas temperature at densities 
above fO^ cm^"^. Future studies can then perhaps consider 
the hydrostatic equilibrium of "Bonnor-Ebert spheres" for 
a more reasonable temperature distribution. 

A brief summary of the results presented here as well 
as a more general discussion of the physical properties of 
pre-protostellar cores is given by Walmsley et al. (2001). 
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Appendix A: The homogeneous spherical cloud 

Consider a spherical, homogeneous cloud of radius R, exposed 
to an interstellar radiation field of mean intensity J^. Let P be 
an internal point at a distance r from the center O, and define 
a coordinate system centered in the center of the sphere with 
the z axis along OP. Let 6 be the angle between the z axis and 
a generic direction in space originating from P. Considering 
only absorption by dust, the mean intensity in P is 



Table B.l. Parameters for the interstellar field 
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where t^, (r, 6) is the optical depth of the cloud at frequency 
f along a path from the point P to the edge of the cloud 
in the direction 6. Since the cloud is homogeneous, Ti,(r,9) is 
proportional to the distance from P to the surface of the cloud, 



T^ix, 9) = r^c(vl — x'^ sin^ — x cos 9), 



where t^c is the optical depth at frequency v to the center 
of cloud, and we have set x = r/R. With the substitution 



77 = yl — a;2 sin^ 6 — xcos6, eq. (A.l) becomes 



Mx) = ^ 

4x 



(1 



rj-^ 



-dri + 



e 'dr] 



range 


X, 


P 


W, 


T, 
(K) 


V-NIR 


0.4 ^m 
0.75 ^m 
1 jj.m 
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1 X 10-" 
1 X 10-^^ 
4 X 10-^^ 


7500 
4000 
3000 


MIR 


100 ^m 


-1.8 


5 X 10"' 




FIR 


140 ^m 


1.65 


2 X 10"* 


23.3 



CBR 1.06 mm 1 



2.728 



Table B.2. Parameters for the opacity 



wavelength range Aa 



(c 



)u at Aa 

Hr 



0.1-10 ^m 
10-400 /im 
400 ^m-10 mm 

and 



1 /im 3.9 X 10^ 

140 ^m 1.5 X 10" 
1.06 mm 3.3 x 10" 



-24 



-26 



1.4 

1.6 
2.0 



J4l) = ^(l + 2r.c-e-2^-). 
4r^c 



(A.6) 



The value of J„ at the edge of the cloud varies between ^J'J 
for T^c -^ 00 and J'J for t^c -^ 0. 

Appendix B: Parameters for the interstellar field 
and the opacity 

We have (for the purpose of deriving the analytic formulae 
given in section 3.2) made approximate fits to the interstellar 
radiation field given by Black (1994) as well as to the grain 
opacities of Ossenkopf & Henning (1994). We approximate the 
radiation field by a sum of modified Planck functions given by 
eq. IM) and the opacity by piecewise power laws as follows. 
For each black body (or modified) co mpo nent of the inter- 

the values of the 



B.l 



stellar radiation field, we give in Table 
parameters Xp, p, Wi, Ti discussed in section 3.2. 

The opacity parameters a and Q^ fo r the different fre- 
quency ranges are also given in Table 



B.2 



These are based on 



the fact that one gets a rough fit (20 percent) to the Ossenkopf 
(A. 2) & Henning (1994) opacity results with a piecewise power law 
fit having an index 1.4 below 10 /im, 1.6 between 10 /im and 
400 /im, and 2 longward of 400 /im (opacity proportional to 
v"). The parameters for dust emission (Qe and a^) used in the 
analytic formulation are appropriate for dust grains at roughly 
10 K and we therefore use the v alues given for the wavelength 
range 10-400 ^m in Table p!2|. 



.(A.3) 



Defining the auxiliary variables x± = Tvc{1 ± x), the integral 
can be expressed in terms of exponential-integral functions, 



Jv{^) = 



IIS 



2{x+ 

X + E2{x-) - 



-E2{x+) +x^Eo{x- 



where Eo{x) = e ^ jx. The values of Jv(x 
at the edge of the cloud are, respectively 

J,.(0) = 3':e-^-% 



x^E^{x^)\ , (A.4) 
at the center and 



(A.5) 
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